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Abstract 

We study the scaling behaviors in the wind velocity time series collected at 
the atmospheric surface layer and compare them with two commonly used 
cascade models, the truncated stable distribution and the log-normal model. 
Results show that although both models can describe the change of prob- 
ability density functions from non-Gaussian to Gaussian like distributions 
with the increase of time scale, they can not fit the scaling behaviors ob- 
served in the probability of return and in the moments at the same time. 
This work provides some clues on the understanding of cascade process in 
the atmospheric turbulence. 

Keywords: Atmospheric turbulence, Cascade process, Scaling behavior, 
Truncated stable distribution, Log-normal model 



1. Introduction 

The probability density functions (PDFs) of the atmospheric turbulent 
velocity increments have been found to go from Gaussian like behavior at a 
time scale of few days to stretched exponential like behavior at a time scale 
of one hour This behavior is similar to that observed in fully developed 
homogeneous and isotropic turbulence, where it has long been recognized that 
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the energy is transferred from large scales to smaller scales by some cascade 
process. Although the atmospheric turbulence is obviously non-homogenous 
and non-isotropic, its statistical similarity with homogeneous and isotropic 
turbulence suggests that there may be also some cascade process associated 
with the energy transfer between synoptic scales and smaller scales in the 
atmospheric turbulence. 

To our knowledge, two models are commonly used to describe the PDFs 
of atmospheric turbulent velocity increments. One is the log-normal model 
which was proposed to account for the intermittency in the local homoge- 
neous and isotropic turbulence [2j. The other is the truncated stable distri- 
bution (sometimes called truncated Levy distribution), which was originally 
proposed to resolve the paradox between infinite variance of stable distri- 
bution and finite variance of real economic systems 0, H| . Some researches 
show that both models are also good enough to fit the PDFs of velocity in 
atmospheric turbulence . However, as we will see in this study, the 

underlying cascade processes of the two models are radically different while 
they are mathematically equivalent for data fitting. 

In addition to PDFs, some other functions can also be used to reflect the 
underlying cascade process from the aspect of scaling. In homogeneous and 
isotropic turbulence, the moments of turbulent velocity increments between 
two points separated by a time or spatial interval (called time or spatial 
scale in this paper) are always used as a tool to reveal the underlying cas- 
cade process. Generally, the gth-order moments vary as a power function of 
scales and their power exponents are concave with respect of q. An another 
function usually discussed in the stable or truncated stable distributions is 
the probability of return, which are also found to be a power of scales when 
scales are small. In fact, this scaling behavior reflects a self-afline cascade 
process in the stable or truncated stable distributions at small scalesjH, 0|. 

The aim of this paper is to analyze the scaling behaviors (i.e., the proba- 
bility of return and the moments as a power of time scales) in the time series 
of atmospheric turbulent velocity increments, by using a large amount of 
experimental data collected in the atmospheric boundary layer (see Sec. [3]). 
Then, the results are compared with above mentioned PDF models each with 
a special cascade process. The inconsistencies between these models and ob- 
served scaling behaviors are discussed in detail (see Sec. HJ). It hopes that 
our research will give some clues about cascade process in the atmospheric 
turbulence. 
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2. PDFs of velocity increments 

Before studying scaling behaviors possibly existing in the atmospheric 
turbulent velocity increments, we first analyze their probability density func- 
tions (PDFs) and compare the results with models. Data used in this study 
is collected at a site located in a steppe of northeast Xilinhaote city, in Inner 
Mongolia, China, where wind velocities are measured at a height of 30 m by 
means of a sonic anemometer with a sampling frequency of 20 Hz (Campbell 
CSAT-3). A velocity time series of continuous 320 hours observations (i.e., 
from 13 October 2009 until 26 October 2009) is chosen to analyze PDFs. 
Thus, as many as 2 x 10 7 samples are used here, which can give a better 
estimate of PDFs especially their tails. 

In the sequel, v(t) will denote the modulus of a wind velocity vector while 
it is used to analyzing the dimensionless value v(t): 

m = vM , (i) 

a v 

where a v is the standard deviation of v(t). The increments of normalized 
wind velocity between two points separated by a time interval of r (i.e, by a 
time scale of r) are then defined by 

Av T = v(t + r)-v(t). (2) 

In fact, the non-overlapping velocity increments can also be estimated by the 
sum of velocity increments at smallest time scale: 

t/tq 

Av T = Y^[v(t + tr ) - v(t + 1)t )]. (3) 

i=l 

where r is the smallest time scale and equals to the reciprocal of sampling 
frequency. 

Standardized PDFs of Av T at different time scales are shown in Fig. [TJ 
One can see that these PDFs are almost symmetrical and go from Gaussian 
like behavior at a time scale of more than ten minutes to some non-Gaussian 
behavior (with longer tails than Gaussian) at smaller time scales. This re- 
sult is similar to that in Ref. [lj (see Fig. 4a in this article) where a longer 
time series (from 1960 to 1999) but with lower sampling frequency (1 hour) 
collected at the atmospheric boundary layer has been used. Moreover, this 
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Fig. 1. Comparison of standardized PDFs of velocity increments Av T at different time 
scales (circles) with (a) the truncated stable distributions (lines) and (b) log-normal PDF 
model (lines). From top to bottom, time scales which range from 0.05 s to 1638.4 s increase 
as a geometric sequence with a common ratio of 2. The symbols <tav t is the standard 
deviation of Av T . Plots have been arbitrarily shifted for illustration. 



deformation behaviors of PDFs in the atmospheric turbulence are also simi- 
lar to those observed in the local homogeneous and isotropic turbulence (see 
Fig. 4b in Ref. [lj). In the latter, it has long been recognized that the energy 
is transferred between various scales by the cascade process. It suggests that 
the energy transfer between atmospheric mesoscales and more smaller scales 
may be also related to some cascade process, although the turbulence at these 
scales is non-homogeneous and isotropic and is influenced by many factors 
like frictional drag, evaporation and transpiration, heat transfer, pollutant 
emission, terrain induced flow modification and so on 

In this study, two commonly used PDF models will be screened. One is 
the truncated stable distribution and the other is the log-normal PDF model. 
The truncated stable distribution was first proposed in Ref. 0] and then was 
used to account for the scaling behavior in the dynamics of an economic 
index Q. After that, many other truncated stable distributions have been 
proposed lCH13j . The distinction between these distributions is mainly focus 
on their tail distributions. In this research, we will focus on the one proposed 
in Ref. 10], which has the fewest tuning parameters and a smooth form at 
the same time. In Refs. 0, 0], this truncated stable distribution has been 
reported to fit the PDFs of atmospheric turbulent velocity. This distribution 
dose not have an analytical expression but has an analytical characteristic 
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function $(Jfe) [10|: 
]n$(Jfe) = 



r 



cos(iTa/2) 



X a - (k 2 + X 2 ) a/2 cos ( a arctan 



when < a < 1, and 
ln$(Jfc) = 



1 — i/3 tan ( a arctan 



7 C 



COs(7TQ!/2) 



A Q - {k 2 + X 2 ) a/2 cos (a arctan 



(4) 



1 — i/3 tan ( a arctan 



-ia$\ a - l k 



(5) 



when 1 < a < 2. Symbols a, /3, 7 and A represent four independent param- 
eters which are normally obtained by fitting Eqs. (j4j) andj5]) with the data. 
Equations (jlj) and §5§ can be written in a simpler form [141 ] : 

ln$(Jfe) = CT(— a)[g(A + i£;) a + p(X — ik) a — X a ], (6) 

when < a < 1, and 

ln$(ife) = CT(-a)[q(X + iA;) a + p{X - ik) a - X a - iaX a ~\q - p)k], (7) 

when 1 < a < 2. The parameters are C = — 7 a /[cos(«7r/2)r(— a)], p = 
(1 + 0)/2 and g = (1 — /3)/2. For standardized and symmetrical PDFs, (3 = 
and 0]: 

CT(-a)aA Q " 2 (a - 1) = 1. (8) 

The log-normal model was originally proposed to account for the PDFs of 
velocity increments in the fully developed homogeneous and isotropic turbu- 
lence [2J. This model does not have an analytical PDF either. Its symmetrical 
PDF is: 



1 



2ttA t 



cxp 



Av 2 
2a 2 



exp 



ln 2 (a/cr r )\ da 



2X1 



a" 



(9) 



where a T and A T are parameters depending on scale r. For standardized 
PDF, there is only one independent parameter and the other parameter can 
be obtained by resolving the equation: 

a 2 exp(2A 2 ) = 1. (10) 
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Figure [T] shows the comparison of experimental PDFs with above two 
models. Results show that the truncated stable distribution seems as "good" 
as the log-normal model when fitting to the data. They both fit to the 
experimental PDFs well except their far tails where the data seem to decay 
faster than the models with unknown reasons (see also Fig. 2 in Ref. [5], 
where a similar deviation from log- normal model has been observed). This 
reason may be physical, but also may be artificial. Due to the limited data, 
tail distributions may be underestimated. In short, by only fitting the PDFs 
we cannot say which model is better to account for the data. We also cannot 
say anything more about the underlying cascade process in the atmospheric 
turbulence. This is the motivation for us to look into the scaling behaviors 
in the atmospheric turbulence. 



3. Scaling behavior 

3.1. Probability of return 

Probability of return is defined as the probability density at the origin 
for symmetrical probability density function (PDF). In the data analysis, a 
probability density defined on a small interval (—h,h) can be used as an 
approximation of the probability of return: 

f T (Av = 0) « f T (-h <Av< h), (11) 

if the threshold h is small enough. In Fig. [2J circles denote the probability of 
return of atmospheric turbulent velocity with a threshold of 0.01. At small 
scales, this function behaves as a power of time scales r: 

f T (Av = 0) OC T~ P . (12) 

A least-square fitting shows that the exponent p ~ 0.2145. For comparison, 
the probability of return of Gaussian distribution is also shown here: 



GV(0) = 




where the standard deviation o~ T = o"A6 T - The difference between / T (Au = 0) 
and G T (0) decreases with the increase of r, which implies a convergence to the 
gaussian distribution at large time scales. We also compute the probability 
of return with a much smaller threshold of 0.001 and find that above results 
do not change any more except that there are some statistical errors at large 
time scales and a little larger slope at small time scales. 
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Fig. 2. Probability of return for the wind velocity increments as a function of time scales 
t (circles for threshold h = 0.01 and squares for h = 0.001). For comparison, the proba- 
bility of return for Gaussian distributions (see Eq. (|1 3[) ) is also shown in this plot (black 
dots). The slope of the best-fit line (fitting data when r < 204.8 s) is -0.2145 and the 
corresponding correlation coefficient is about 0.999. Broken line shows a power function 
with an exponent of -0.5. 



3.2. Moments 

Scaling behavior in the statistical moments of velocity increments (also 
called structure functions) is the central topic for the homogeneous and 
isotropic turbulence. Many theoretical works have devoted to the under- 
standing of this scaling behavior and the underlying physical mechanisms 



15| . For the non-homogeneous and non-isotropic atmospheric turbulence, 
its statistical moments also have scaling behavior at a wide range of time 
scales. Figure [3^ shows the statistical moments of atmospheric turbulent 
velocity increments. One can see that at small time scales, 



Z T) g = \Av T \i oct*«, (14) 

where the exponents £ g vary as a function of q (Fig. |3b). This power-law 
scaling behavior is observed for time intervals spanning three orders of mag- 
nitude for small values of q. However, for larger values of q (for example, 
q = 2.5 and 3 in Fig. |3^) data begin to deviate from this scaling behavior at 
very small time scales (r < Is). The higher-order moments is mainly defined 
by the tail distributions which describe the statistical behavior of turbulent 
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Fig. 3. (a) Scaling behavior of qth order moments Z Ttq . Lines represent the least-square 
fitting to Eq. (fT4|) and the corresponding correlation coefficients R are also listed. Grey 
points are not used in the fittings, (b) Scaling exponents £ q vary as a function of q. Line 
shows the fitting result with the log-normal model and the broken one is a line through 
the origin with a slope of p, where p is the scaling exponent of probability of return (see 
Eq. (HID). 



eddies with large characteristic velocities. It suggests that the relatively very 
small turbulent eddies with large characteristic velocities may have different 
cascade process from these eddies with larger scales. However, as already 
remarked in Fig. (TJ all the conclusions about the tail distributions should be 
carefully screened because of the limited data. 

Comparing with the exponents £ 9 of homogeneous and isotropic turbu- 



lence [15], we found that £ 9 of non-homogeneous and non-isotropic atmo- 
spheric turbulence are also concave but their values are much smaller. For 
the homogeneous and isotropic turbulence at high Reynolds number, the 
slope of £q at small values of q is about 1/3 and £3 = 1 (Kolmogorov's four- 
fifths law). As shown in Fig. 12b, the slope of £ q at small values of q in 
the atmospheric turbulence is about 1/5 and £3 ~ 0.35. According to the 
Kolmogorov hypothesis in the log-normal model^, fl6| . 

9A 2 = Aq - yU In r, (15) 

where Ao and /1 are parameters independent of r. For a fixed time scale, 
a smaller value of fi corresponds to a larger value of A which means that 
the tail of PDF becomes longer and the corresponding turbulence is more 
intermittent (see Eq. Following Eq. ([15]) and the assumption o T oc r b , 
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one can deduce that the exponents £ g behave as a simple parabola: 

i q = bq- aq 2 , (16) 



where b is some constant and a = fi/18. A least-square fitting shows that a ~ 
0.04 for atmospheric turbulence while this value is about 0.01 observed in the 
fully developed homogeneous and isotropic turbulence [l5j. Thus, from the 
view of log-normal model and the Kolmogorov hypothesis the homogeneous 
and isotropic turbulence is less intermittent than the non-homogeneous and 
non-isotropic atmospheric turbulence. 



4. Comparison with models 

The truncated stable distribution and the log-normal model have different 
cascade processes which can cause the different scaling behaviors. The center 
part of truncated stable distribution at small scales behaves like a stable 
distribution with the same parameters a, (3 and 7. A notable feature of the 
stable distribution is its stability [171 ]: 



Ati r0i! = n*Av T0 , (17) 

i=i 

where n is a positive integer and A-0 rOii (i = 1,2, ■■■,n) are independent 

random variables with the same distribution as Av TQ . The symbol "=" means 
equality in distribution. 

From Eq. fl3]), one can see that 

n 

Av T = A<i, (18) 



i=i 



where Av TQti = v(t + ir ) —v(t+ {i — l)r ) and n = t/t . Thus, the truncated 
stable distribution contains a self-affine cascade process at small scales jsj: 

Av T = niAv TQ . (19) 
An equivalent expression of Eq. ( |T9l) is 

f n) (Av)=n 1 / a f T (n 1 / a Av). (20) 
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From this expression, one can deduce that the probability of return for the 
truncated stable distribution behaves as a power of time scale: 



f T (Av = 0) OC T' l/a . (21) 

Although the probability of return of truncated stable distribution has a 
scaling behavior at small scales, one should note that the scaling exponent 
1/a G [0.5,oo). This scaling exponent is larger than that observed in the 
atmospheric turbulence (see Fig. [2]). For the gth-order moments, authors 



(Refs. [14|, |l8j) have proved that the truncated stable distribution also has a 



scaling behavior but the scaling exponents vary as a bi-linear function of q: 

f q/a, (0 < q < a) , . 

* ~ I 1, (? > «) { } 

which is significantly different from the scaling exponents observed in the 
atmospheric turbulence, as shown in Fig. [3b. 

In order to deduce Eq. flUJ), we have assumed that wind velocity incre- 
ments in Eq. (|T7j) are independent of each other . However, some correlations 
may exist among these increments. A direct proof comes from the scaling 
behavior of second-order moment. As shown in Fig. [3j the scaling exponent 
of second-order moment is obviously smaller than one, while for independent 
random variables this value equals to one. Correlation may be a possible 
reason that the scaling exponents of probability of return is smaller than 
the prediction of truncated stable distribution. To verify this speculation, 
we generate an independent time series by randomly sort the original time 
series of Av T0 . This surrogate time series will have the same statistical dis- 
tribution as the original time series. Independent surrogate time series at 
other larger scales can also be obtained by Eq. ([3]). We then compute the 
probability density functions of the surrogate time series at different time 
scales and find that with the increase of scales they approach Gaussian-like 
distributions much faster than the original time series (comparing Fig. 
with Fig. [TJ. This is also proved in Fig. Hb, where the probability of return 
is almost the same as that of independent Gaussian distribution when time 
scale r > Is. Most importantly, one can see that the scaling exponent of 
probability of return is about 0.78 which falls on the range [0.5, 00) that 
the scaling exponent of truncated stable distribution belongs to. Thus, we 
can conclude that the correlation in the wind velocity increments is indeed 
a possible reason that the scaling exponent of probability of return deviates 
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Fig. 4. Results for the surrogate wind velocity increments without correlations: (a) Stan- 
dardized PDFs at different time scales. From top to bottom, the corresponding time scales 
are the same as those in Fig. [TJ Plots are also arbitrarily shifted for illustration, (b) Prob- 
ability of return with threshold h = 0.01 (circles). The probability of return for Gaussian 
distribution (see Eq. (1131) ) is also shown here (black dots). The slope of the best-fit line 
(fitting data when t < 0.8 s) is —0.7808 and the corresponding correlation coefficient is 
about 0.99. Broken line shows a power function with an exponent of -0.5. (c) Scaling 
behavior of g-th order moments Z T>g . All the symbols are the same as those in Fig. [3^. 
(d) Scaling exponents £ g vary as a function of q. Broken lines shows a bi-linear function. 
The slope of this function is 0.7808 when q < 1/0.7808 (see Eq. ([22])). 
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from the prediction of truncated stable distribution. However, this reason 
is not enough to account for the deviation of moments. For the surrogate 
time series, their exponents of moments £ g is also a concave function of q, 
while the scaling exponents of truncated stable distribution obey a bi-linear 
behavior. 

According to Eq. the log- normal PDF model can be expressed by a 
simple random mapping: 

AfV = W 9r)Xr G, (23) 

where Wa T ,x T is a log- normal random variable and its logarithm has a mean 
of lnov and a standard deviation of A T . The symbol G denotes a normally 
distributed random variable and is independent of W& Tt \ T . For log- normal 
random variables, we have 

W- aTiXT = W- CTiXT W- aThKn (24) 

where Wa Tl \ T and W„ , are independent log-normal random variables, and 
In or + hi ay = lna T and \\ + \ 2 T , = X 2 .. Based on Eq. ([23]) and Eq. f[2ij) . 
one can deduce that the log-normal PDF model has a cascade process: 

Av T = W 5tM /\v t , (25) 

Comparing Eqs. ( TT§j) and (12"51) . we find that the cascade process is very differ- 
ent between the truncated stable distribution and the log-normal PDF model. 
For the truncated stable distribution, the connection between different scales 
is a non-random power function which can only produce a stochastic process 
with scaling exponents £ q of moments varying as a linear function of order 
q. This process is referred as a self-affine or self-similar fractal. To obtain 
a stochastic process with a non-linear scaling exponents of £ q which is re- 
ferred as a multifractal, one should extend the non-random power function 
to a random variable [8| . For the log-normal PDF model, the connection is a 
log-normally distributed random variable and it can produce a multifractal 
process by assuming suitable relationships between parameters (a r , A r ) and 
time scales r. 

As we have mentioned in Sec. I3.2[ under the assumptions of Eq. f lT5|) and 
a T oc r b , the scaling exponents £ g of log- normal model vary as a parabola 
function (see Eq. f lT6|) ) and this function can fit the data well (see Fig. |3b). 
However, we will show that whatever the relationships between parameters 
(o" T , A T ) and time scales are assumed the log-normal PDF model can not fit 



12 



1.6 



< 



1.5- 
1.4 - 
1.3 - 
1.2 - 



0.9 - 
0.8 



0.74 
10" 



10"' 10" 



10 

r(s) 



10' 10" 



10" 



Fig. 5. Test of Eq. (|29j) . 

the probability of return and the moments at the same time. According to 
Eq. (jUJ), the probability of return of log-normal PDF model is: 



fr(Aij = 0) 



2na T 



We can also deduce its g-order moment [2J: 

Z r , q = a«P q exp{q 2 \ T 2 /2), 

where 

1 r°° 

P q = —= / |x| 9 exp(-x 2 /2) dx. 

V27T J-oo 

From Eqs. f )26|) and fT2T|) . we have: 

7T 2 ^ |1 / r (Al) = 0) i 

2 Z r ,2 



(26) 

(27) 
(28) 

(29) 



In Fig. [51 one can see that the left-hand side of Eq. (1291) is not a constant 
of one and varies as a function of r. Thus, we conclude that the log-normal 
PDF model can not fit the probability of return and the moments at the 
same time. 
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5. Conclusions 



In this paper, we have shown the scaling behavior of the probability 
of return and the moments of atmospheric turbulent velocity increments. 
The probability of return is found to vary as a power function of time scale 
at smaller time scales. This scaling behavior is observed for time scales 
spanning three orders of magnitude, from 0.1s to 100s. The moments are 
also found to vary as a power function of time scales at smaller time scales. 
The scaling range can also spanning three orders of magnitude when the 
order of moment is not too large. However, with the increase of the order 
of moment the scaling range begins to shrink. Since the behavior of higher- 
order moments is mainly defined by tail distributions, the shrinkage may be 
caused by different cascade processes that the central and tail distributions 
represent. This conclusion should be carefully screened by using more data. 
Scaling exponents of moments vary as a concave function of order, which 
suggests that the time series of atmospheric turbulent velocity increments is 
a multifractal. Scaling exponents of moments observed in the homogeneous 
and isotropic turbulence are also concave but the values are obviously larger 
than those observed in the atmospheric turbulence. 

The scaling behavior of the probability of return and the moments ob- 
served in the time series of atmospheric turbulent velocity increments is com- 
pared with the predictions of two commonly used PDF model, the truncated 
stable distribution and the log-normal PDF model. We find that both PDF 
models can not fit the probability of return and the moments at the same 
time. The truncated stable distribution can produce a power law of proba- 
bility of return but the power exponent is too large to fitting the data. One 
possible reason is that the truncated stable distribution can only describe 
the independent random process while some statistical correlations always 
exist in the time series of atmospheric turbulent velocity increments. The 
truncated stable distribution can also produce a power law of moments. The 
power exponents vary as a bi-linear function of order, while the observed 
exponents are concave. The log-normal PDF model seems to be better than 
the truncated stable distribution. Under some conditions, its exponents vary 
as a parabola function which can fit the data well. However, this model can 
not fit the probability of return and the moments at the same time either. 

Although above mentioned PDF models can not describe the scaling be- 
havior in the atmospheric turbulence, they also give a clues that the scaling 
behavior may related to the cascade process of turbulent eddies. Our study 
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may help in the understanding of the cascade mechanism in the atmospheric 
turbulence. 
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